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ABSTRACT 

India’s first Moon mission Chandrayaan-1 carrying eleven scientific instruments was launched on 22nd 
October 2008. Lunar Laser Ranging Instrument (LLRI) is one of the scientific instruments carried by 
Chandrayaan-1. It is an instrument aimed to enhance the study of the Moon’s surface. LLRI instrument 
was commissioned on 16 th November 2008. The ISRO’s mission operational orbit determination (OD) 
software was suitably updated to estimate the Chandrayaan-1 orbit with LLRI measurements. This 
paper describes the method involved in OD process and orbit solutions with simulated as well as live 
LLRI data of Chandrayaan-1. This work resulted in obtaining satisfactory results with LLRI 
measurements. 

1. INTRODUCTION 

India’s first Moon mission Chandrayaan-1 carrying eleven scientific instruments for the purpose of 
expanding scientific knowledge about the Moon was launched on 22 nd October 2008 from Satish 
Dhawan Space Centre, Sriharikota, India by India’s Polar Satellite Launch Vehicle PSLV-C11. 
Spacecraft was injected into transfer orbit of (254.4 X 22932.7) km with inclination of 17.9 deg at 
2008-10-22-01-10-19-081 UT. The main objective of the mission is a simultaneous chemical, 
mineralogical and photo geologic mapping of the whole Moon with high spatial resolution using high 
resolution state of the art sensors. 

The spacecraft was put into Moon’s polar, circular orbit of about (100 X100) km by carrying out five 
Earth bound maneuvers (EBNs), one trajectory correction (TCM) and four Lunar bound maneuvers 
(LBNs) on 12 th November 2008. Lunar Laser Ranging Instrument (LLRI) is one of the eleven 
scientific instruments carried by Chandrayaan-1. It is a pulsed laser-ranging instruments aimed at 
enhanced study of the Moon’s surface. LLRI measurements can be obtained from both the dark and 
sunlit portion of the Moon, thereby significantly increasing the useful observational coverage. LLRI 
instrument for Chandrayaan-1 capable of making lunar topography measurements with a resolution of 
less than 5 meters was successfully turned on 16 th November 2008 at 03:47:088 UT. The purpose of 
this paper is to study the feasibility of Chandrayaan-1 orbit determination using LLRI measurements. 

ISRO’s mission operational orbit determination (OD) software [1] was suitably updated to determine 
the orbit with LLRI measurements and demonstrated using simulated as well as actual LLRI 
measurements obtained during initial phase of the mission. 

The main features of orbit detennination are Trajectory Generation, Observation Modeling and 
Estimation. The main difference between OD with tracking data and with LLRI data is observation 
modeling. This paper describes the method involved in observation modeling in detail and OD results 
with simulated as well as actual/live LLRI data. The comparison study was made between OD results 
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using tracking data (range and Doppler data) and using LLRI data. The LLRI based OD results will 
substitute in case of either non-availability of tracking data or non-availability of converged state with 
tracking data. The work on updating the Chandrayaan-1 operational orbit determination program using 
tracking data for LLRI measurements resulted in very satisfactory OD results. 

2. ORBIT DETERMINATION SYSTEM 

A precise knowledge of orbital parameters is a prerequisite for the determination of the position of 
satellite at any given time. Mathematically orbit is specified by a set of six orbital parameters which 
completely describes the motion of the satellite within the specified accuracy over a period of time. 
Determination of the satellite orbit means the refinement of the initial orbit parameters. The problem of 
OD consists of comparing measurements taken on a satellite trajectory to a model representing that 
trajectory. The model is generally represented by a system of differential equations whose constants of 
integration define the satellite trajectory as a function of time. Thus, given a priori estimate of these 
constants of integration (state parameters), the problem is to update or correct the initial (a priori) state 
parameters as a function of measurements. The orbit detennination system consists of tracking data 
system and data handling, dynamic models, computational techniques, which include trajectory 
generation and estimation technique. 

The main computation process of the orbit determination system is trajectory generation, observation 
modeling and estimation. Trajectory generation is performed through numerical integration of 
differential equations of motion of satellite. The force model for artificial satellites normally includes 
gravitation attraction of all planets, aerodynamic drag, luni-solar gravitation and solar radiation 
pressure. Cowell's method is used for trajectory generation. Weighted least squares technique and 
iterative differential correction process is used to obtain the refined state. The method of integration of 
equations of motion and estimation are described in the following sections. 

2.1 Main features of Orbit Determination System 

The main features of orbit determination system include measurement data processing, trajectory generation, 
observation modeling and orbit estimation. Generally observational data is tracking data. The present study of 
obtaining orbit determination is based of LLRI measurements. 

2.1.1 Method of Integration for Trajectory Generation 

The coupled nonlinear second order differential equations of motion are integrated numerically through 
Cowell's method. The chief advantage of special perturbation technique is its high accuracy. The 
perturbative acceleration acting on the satellite is modeled. Numerical integration methods are 
customarily classified into two parts namely single step and multi step methods. Multi step methods 
normally need a starter. The single step methods that are considered are those of Runge Kutta or RK- 
Gill etc. Multi step methods are "Adams" or its variants. Exclusively for solving second order 
equations, double integration method is the recommended procedure. These methods are stable. Gauss- 
Jackson-Merson’s (GJM) 8 th order method is employed here [2]. 

The basic equations to solve the equation of motion 


x = f(x,x,t) ( 1 ) 

are 
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v = x(t. + h) = h 2 (V - 2 Xi+—Xi+—Vx t +—V 2 Xi+— V 3 x,+V 4 x,+...) 

' 12 12 240 40 12096 (2) 

1 5 ^ 951 95 190R7 

x M =x(t i +h) = h(V~ 1 x+-Xi+—Vxi+-V 2 Xi+ -V 3 x,-+-V 4 x,+-V 5 x,+...) 

2 12 8 720 288 60480 


in standard backward difference notation. 


The Gauss-Jackson corrector fonnulae are 


= /? 2 (V 2 X; + — Xi+i-— V 2 Xi+i -— V 3 Xi+i- 

12 240 240 60480 
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xm =h(V 1 Xi+ — x, +1 - — Vxi+i—— V 2 x, +1 —— V 3 x ; +i—^-V 4 x,+i—V 6 Xi+i)+.„.) (3) 

2 12 24 720 160 24192 

2.1.2 Observation Modeling 

For a given LLRI measurement, the theoretical/computed measurement p c is obtained using satellite 
definitive ephemeris and spacecraft quaternions (q s ) along with instrument misalignment angles. 

The following steps are followed to obtain theoretical measurement: 



The observed LLRI measurement p 0 , the scalar value is converted into inertial vector p 0 , using direction 
cosine (DC) matrix obtained from spacecraft quaternions and instrument misalignment matrix. 

p 0 = p 0 [ DC matrix obtained from q s ] [instrument misalignment matrix] 


The vector R, from Moon’s centre to Moon’s surface is given as R = r + p 0 , where r is the spacecraft 
position vector from definitive ephemeris (inertial EME-J2000). R is converted from selenocentric to 
selenographic to compute longitude X and latitude (p of LLRI measurement on the lunar surface. 

The Moon’s radius R at a given longitude and latitude is calculated using the following relation 

R=XX Pi m (sin(p)(C| m cos(mk)+ S| m sin(im)) (4) 

Where, Pi m is the Lagrangian coefficient of order 1, m; 

(p and X are the latitude and longitude; 

Cim, Si m are topographic coefficients 


The computed R is converted into vector R using X and (p and converted into inertial (EME J2000) Rj. 
Computed LLRI measurement p c using Rj and r is p c = Rj - r . 
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2.1.3 Estimation 


Optimal estimate of satellite state is computed from a given vector of tracking measurements which 
depends on the instantaneous position and velocity of spacecraft. The estimated state is valid over the 
period during which measurement data is collected. 

The general procedure for all definitive orbit computations is to set up some dynamical model of the 
orbit and use the observations to improve the orbit parameters of the model by the process of 
differential correction. Weighted Least Squares estimation process is applied. Observations are 
considered for some pre-selected time period and by differential correction of parameters of the model, 
the sum of the squares of the residual is minimized. The basic idea of least-squares estimation as 
applied to orbit detennination is to find trajectory and the model parameters for which the square of the 
difference between the modeled observations and the actual measurements becomes as small as 
possible. In this estimation process of “Weighted Least Squares” it is necessary to compute partial 
derivatives of observations with respect to model parameters. Brief description of method of estimation 
adopted in this present study is as follows. 

Given m observations and an orbit that is to be corrected using these observations consisting of n 
parameters, Pj,j = \,...n, and if M o is the observed quantity, M c is the corresponding computed 

quantity from the mathematical model, then each of the observed quantity can be considered as a 
function of state parameters, the station coordinates, time and if necessary of other relevant parameters 
affecting the motion of the satellite. In this context, only the dependence of orbital parameters is of 
importance. 

Thus, 

M c =M c {P} (5) 


By comparing the computed quantities with m observed quantities, the residuals are given by, 


= M oi -M ci {P.},i = 


( 6 ) 


The object is to obtain expressions in terms of known quantities. The residuals after the first iteration 
become: 


AM' =M oi -M ci {Pj'},i = \,2,...m (7) 

where, 


p j = p j 


A P, 


Expanding by Taylor series and linearizing, the conditional equations become: 
AM',- = AM,. 


JM . 

E- ,j = 1,2, ..m 


J =1 dP, 


( 8 ) 


By multiplying weighting factors, for all m observed quantities the equations of condition can be 
written in matrix notation as: 
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W t R =W l R-W l AAP 


(9) 


where, 

W l is the weighting matrix of mean measurement error (total expected error for the each measurement 
type due to both random noise and systematic errors) which is derived from the covariance matrix of 
observations by taking square root of the diagonal elements. R' is the column matrix of residuals after 
differential correction, R is the column matrix of residuals before differential correction, A is the matrix 
of partial derivatives. 

If everything is ideal, R' should contain all zeros after just one iteration. The solution for A P would 
have been trivial if there were same observed quantities as that of parameters to be refined. In practice, 
'm' is much greater than ‘n’ and the second order terms are not negligible, so that one iteration does not 
suffice. Therefore, weighted least square technique is adopted. The condition of the least square 
approach is: 

\W X R ] T [W X R ] is minimum (10) 

that leads to the nonnal equations: 


A T W 1 r (W 1 AAP-W 1 R) = 0 


(ID 


Letting 

w = w 1 r w l , 

the equation yields: 

AP = (A T WAY l A T WR (12) 


A P is the correction to P for the current iteration. This process is to be carried out for a number of iterations till 
the solution converges. 

3. MEASUREMENTS USED FOR ORBIT DETERMINATION 

Orbit detennination using LLRI measurements is demonstrated with Chandrayaan-1 simulated as well 
as actual/live data. About two days of LLRI measurements were simulated using Chandrayaan-1 state 
parameters given in Table 1. 

Table 1 Measurements Simulation 


Epoch : 2008-12-18-15-48-54-816 

Start of measurement : 2008-12-18-15-49-54-816 
End of measurement : 2008-12-20-15-18-54-816 


State parameter 


(km, km/s, EME-J2000) 


X.Y.Z 

68.485, -278.061. 1828.229 

XD,YD,ZD 

-0.111485 ,-1.603584 , -0.230156 
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Live LLRI data obtained on 16 th November 2008 (1 st day of commissioning LLRI instrument) and 1 st 
week of February 2009 was used to demonstrate OD. Orbit detennination results using 16 th November 
2008 LLRI data are shown in Table 3. These OD results are compared with orbit solutions obtained 
using tracking data. Table 4 gives summary of OD results comparison using 1 st week of February 2009 
LLRI data. 

4. OD RESULTS AND COMPARISON 

Orbit detenninations were carried out with simulated and actual Chandrayaan-1 LLRI data and 
comparison study was carried out with OD results obtained using tracking data (Range and 
accumulated Doppler data) 


Table 2 Orbit determination summary using simulated data 


EPOCH 


: 2008 12 18 15 48 54 816 


TRACKING DATA FROM 
TO 

TRACKING DATA DURATION 
No. of LLRI Measurements 


: 2008 12 18 15 49 54 816 
: 2008 12 20 15 18 54 816 
: 47.48 hours 
: 2850 


Definitive State Parameters (EME-J2000) 


Converged Versus Actual State (km. km/s) 


Actual : 

Realized: 

Difference (Actual - Realized) 


68.485 -278.061 
68.471 -278.111 
0.014 0.05 


1828.229 -0.111485 
1828.222 -0.111467 
0.007 -0.000018 


-1.603584 -0.230156 
-1.603579 -0.230201 
-0.000005 0.000045 


Table 3 Orbit determination summary using actual LLRI Data 


EPOCH 

TRACKING DATA FROM 
TO 

TRACKING DATA DURATION 


: 2008 11 16 00 02 00 000 
: 2008 11 16 04 28 48 401 
: 2008 11 16 17 23 39 442 
: 12.91 hours 


Definitive Orbit (EME J2000) 


OD with tracking data. 
OD with LLRI data : 
Difference : 


Orbit Comparison 

a (met) e i (deg) to (deg) Q (deg) 

1842015.732578 0.005425 93.106987 260.379444 86.307290 
1841976.110290 0.005396 93.342136 261.742166 86.517621 
39.622288 0.000029 0.235149 1.362722 0.210331 


M(deg) 
304.758661 
303.472880 
1.285781 
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Table 4 Orbit comparison (LLRI OD Vs. Tracking data OD) 


Epoch 

Orbit with Trk.data 
‘a’(km) V ‘i’(deg) 

Orbit with LLRI data 
‘a’ (km) V ‘i’(deg) 

Difference 

‘a’(km) ‘e’ ‘i’(deg) 

2009 02 01 

1839.94328 0.004194 92.8556 

1839.95064 0.004225 92.9648 

0.00736 0.31e-4 0.1092 

2009 02 02 

1840.65106 0.003560 92.8873 

1840.69025 0.003644 92.8527 

0.03919 0.85e-4 0.0346 

2009 02 03 

1839.43615 0.002965 93.1628 

1839.42036 0.003023 93.2679 

0.01579 0.58e-4 0.1051 

2009 02 04 

1839.94554 0.002327 93.4794 

1839.99136 0.002318 93.4702 

0.04582 0.90e-5 0.0092 

2009 02 05 

1840.33579 0.003244 93.5257 

1840.34211 0.003298 93.4699 

0.00633 0.55e-4 0.0558 

2009 02 06 

1839.28389 0.003485 93.4894 

1839.30033 0.003523 93.2542 

0.01644 0.38e-4 0.2352 

2009 02 07 

1840.41730 0.005942 93.6716 

1840.42967 0.006139 93.6867 

0.01237 0.20e-3 0.0151 

2009 02 08 

1839.46514 0.008093 93.6718 

1839.44635 0.008272 93.7944 

0.01878 0.18e-3 0.1226 


5. CONCLUSION 

The LLRI based OD results can substitute in case of either non-availability of tracking data or non¬ 
availability of converged state with tracking data. The work on updating the Chandrayaan-1 operational 
orbit determination program using tracking data (range and Doppler data) for LLRI measurements 
resulted in very satisfactory OD results. 
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